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actions. A formation zone intranuclear cascade of low energy secondaries inside the target 
nucleus is discussed. We calculate excitation energies of residual nuclei left after the intranu- 
clear cascade process and treat their further disintegration by introducing models for the 
evaporation of protons, neutrons, and light fragments, high energy fission, and by applying a 
Fermi Break-up model to light nuclear fragments. The results are compared to data on target 
associated particle production. We furthermore calculate cross sections for the production of 
nuclear fragments. 
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1 Introduction 



The Dual Parton Model (DPM) |l| and Monte Carlo (MC) implementations of this model 
for hadron-hadron ||^, ^ , hadron-nucleus, and nucleus- nucleus collisions |^ have been quite 
successful in describing many aspects of hadron production in high energy collisions. So far 
however, MC models for hadron-nucleus and nucleus-nucleus collisions based on the DPM did 
mainly describe the high energy component of newly created hadrons, not the many particles 
resulting from the nuclear disintegration following the high energy collisions. 

Models for nuclear evaporation and fragmentation and for high energy fission are however 
usually included in hadron cascade models such as FLUKA ||^, || used for detector sim- 
ulation and for the evaluation of radiation damage to high energy accelerator and detector 
components. In the older models, the nuclear excitation energy, which is the starting point for 
the calculation of the nuclear disintegration, was often introduced only in a phenomenological 



way or it was calculated on the basis of valid but simple intranuclear cascade models ||10|, ^ 
which are not applicable in the multi-GeV energy range of present experiments. The intranu- 
clear cascade models have been greatly improved since then p, |^, ^ and their range of 
validity can be extended to higher energies due to the introduction of the formation zone 
concept [0, [1^ . 



Here, we use the DTUNUC |||] and DPMJET-II g MC implementations of the DPM for high 
energy hadron-nucleus and nucleus-nucleus collisions. These MC models contain a formation 
zone intranuclear cascade which is responsible for knocking out cascade protons and neutrons 
of the residual nucleus. The cascade protons have energies which are typical for the so called 
grey prongs observed in emulsion experiments. Therefore, we are able to calculate the nuclear 
excitation energy of the residual nucleus. In a second step, this excitation energy is the 
basis for nuclear evaporation, and high energy fission reactions. In the present paper these 
mechanisms are investigated in hadron-nucleus collisions. In a forthcoming publication we 
intend to extend these studies to peripheral high energy nucleus-nucleus collisions. 



In Sec. 2 we describe the formation zone intranuclear cascade model ||T^ and focus on the 
calculation of excitation energies. In Sec. 3 models for evaporation and fragmentation are 
presented. In Sec. 4 we compare computed cross sections and multiplicities of grey and black 
prong production to experimental data. Furthermore, cross sections for the production of 
residual nuclei are discussed. In Sec. 5 we summarize our results. 

2 The calculation of excitation energies in the forma- 
tion zone intranuclear cascade model 

2.1 The two-component Dual Parton Model for hadron-nucleus 
and nucleus-nucleus collisions 

The two-component DPM and its MC realizations have been discussed in detail in ^ |^, 
^, ^, |16|. Therefore, we summarize only briefly the main steps leading to the multiparticle 



state, which is the starting point for the intranuclear cascade and evaporation models being 
described in this paper. 
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The MC model for hadron-nucleus and nucleus-nucleus interactions starts from an impulse 
approximation for the nucleons of the interacting nuclei. The spatial initial configuration, i.e. 
the positions of the nucleons in space-time in the rest system of the corresponding nucleus, is 
sampled from standard density distributions. For energies above 3-5 GeV/nucleon the collision 
proceeds via v elementary interactions between and vt nucleons from the projectile and 
target, resp. The values Up, and Ut are sampled according to Glauber's multiple scattering 
formalism using the MC algorithm of [0. The particle production is well described by the 



two-component DPM which is applied as in hadron-hadron interactions 0, IC, 18|. As a 
result a system of chains connecting partons of the nucleons involved in the scattering process 
is formed. The chains are hadronized applying the model JETSET ||19| , pO| . The hadrons 
may than cause intranuclear cascade processes, which are treated by the formation zone 
intranuclear cascade model |]T5|, an extension of the intranuclear cascade model ||T0|, [Tl|. At 
energies below 3-5 GeV/nucleon the formation zone intranuclear cascade model itself provides 
a reasonable description of inelastic nuclear collisions. 

In the following we summarize the main ideas of the formation zone intranuclear cascade 
model for hadron-nucleus interactions. Modifications which have to be introduced to describe 
nucleus-nucleus collisions will be discussed in a forthcoming paper. The physical picture 
explaining the absence of the intranuclear cascade at high energies is the concept of the 



formation zone |T3|. It has been introduced in analogy to the Landau- Pomeranchuk 



effect, which explains the observation that electrons passing through material become more 
penetrating at high energies. For the formation zone of an electron with 4-momentum p and 
energy E upon radiation of a photon with 4-momentum k one obtains 

E E 1 

^ = 1 = , (1) 

k ■ p muje 



where Ue is the frequency of the photon in the rest frame of the electron and E /mis the time 
dilatation factor from the electron rest frame to the laboratory. Within the quark model, the 
states being formed in the primary nucleon-nucleon interaction can be understood as consisting 
of valence quarks only, i.e without the full system of sea quarks, antiquarks, and gluons and 



have therefore a reduced probability for hadronic interactions inside the nucleus [15|. The 



formation zone concept can be translated to hadron production as follows |T^. Denoting the 



4-momenta of the projectile hadron Pp and the secondary hadron pg in the laboratory frame 
with 



Pp = {Ep, 0, 0, ^Ej - ml), p, = {E,,p,^, ^El - ml - f^^) (2) 

and replacing in Eq. ([1|) the electron momentum by Pp and the photon momentum by p^, the 
hadron formation zone reads for Ep ^ 

lEg Eg 
{mpxy + ml + p'^gj_' Ep 

Since for most of the produced secondaries the term {mpxY can be neglected one can approx- 
imate 

^Lab ~ IsTs, Is = —■ (4) 
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In the rest system of the secondary hadron s, we define an average formation time Ts needed 
to create a complete hadronic state |14 



= To 2^2 - (5) 

mj + pi^ 

To is a free parameter, which has to be determined by comparing particle production within 
the model to experimental data. Typical values are in the range from 1 fm/c to 10 fm/cQ. 
From the comparisons discussed further below we find tq = 2 fm/c. For each secondary we 
sample a formation time r from an exponential distribution with an average value as 
given in Eq. (^. As it was described in in our MC model we know the full space-time 
history of the collision. In any particular Lorentz frame we can follow the trajectories of the 
secondaries created in the hadronization of the chains in space and time. Due to relativistic 
time dilatation secondaries with high energies in the nucleus rest system are mostly formed 
outside the nucleus and are therefore not able to initiate intranuclear cascade processes. On 
the other hand, the lower the energy of the secondary hadronic system the higher is the 
probability to form a hadron inside the nucleus. These hadrons may therefore reinteract 
with spectator nucleons taking into account the nuclear geometry. In the model, inelastic 
secondary interactions of hadrons with energies below 9 GeV are described with the code 
HADRIN |2^. In general, the intranuclear cascade would start with resonances resulting from 



the chain hadronization procedure, but we apply the following way. Since the interaction 
cross sections of resonances needed within HADRIN are less well known, we firstly treat their 
decay and sample the intranuclear cascade starting from the stable particles. We assume 
that the effect of this approximation can mostly be compensated by using an effective tq 
value. Reinteractions within the colliding nucleus beyond 9 GeV are very rare and therefore 
neglected in the present approach. Pauli's principle is taken into account as described in 
For the secondaries produced in intranuclear cascade processes we apply the same formalism as 
described above and therefore obtain a formation zone intranuclear cascade in all generations 
of secondaries. 



2.2 The calculation of nuclear excitation energies 

The treatment of nuclear effects within the MC model has already been discussed in . Since 
they are essential in calculating excitation energies of nuclei left after primary interactions and 
intranuclear cascade processes we summarize the basic ideas. Fermi momenta for nucleons 
as well as a simplified treatment of the nuclear potential are applied to control the gener- 
ation of low-energy particles. Nucleon momenta are sampled from zero-temperature Fermi 
distributions 

= AT^^p f . (6) 
dp (Pf^)^ 



^In Refs. 1^, ||] To was fixed to tq 5 fm/c whereas in Ref. jTs] to=1-2 fm/c was used. 
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Here and in the following the indices "n" and "p" denote neutrons and protons, resp. The 
maximum allowed Fermi momenta of neutrons and protons are 



n,p 

Pf 





3h^' 


A Va ) 


87_ 



(7) 



with Va being the volume of the corresponding nucleus with an approximate nuclear radius 
Ra = ''"o^^''^; ''"o = 1-29 fm. Modifications of the actual nucleon momentum distribution, as 
they would arise, for instance, taking the reduced density and momenta in the nuclear skin 
into consideration, effectively result in a reduction of the Fermi momenta as compared to those 
sampled from Eq. (0). This effect can be estimated by a correction factor a^od which modifies 
the Fermi-momenta. Results presented in this paper have been obtained with a^^^=0.75. The 
depth of the nuclear potential is assumed to be the Fermi energy and the binding energy for 
outer shell nucleons 

^''''n,p 

To extend the applicability of the model to the energy region well below 1 GeV an approximate 
treatment of the Coulomb-potential is provided. The Coulomb-barrier modifying the nuclear 
potential is calculated from 

Z1Z2 

" 47reoro (A'/^ + A^^') 

with the mass numbers Ai, Ai and charges Zi, of the colliding nuclei, i.e. with = |Zi| = 1 
for charged hadrons entering or leaving the target nucleus, e denotes the elementary charge 
and ro = 1.29 fm. 

The excitation energy V of the residual nucleus with mass number Ares and charge Zres, 
i.e. the energy above the ground state mass -Eo,res, is given as 

^ -^res -^0,res; 

-E^O,res -^rcs^p ~l~ (^rcs -^rcs 

)mn — i?bmd(^res, -^res)- (10) 

We calculate the binding energy £'bind(Ares, ^res) using the experimentally determined excess 
masses of all known (measured) nuclides and using mass formulae for nuclides far from the 
stable region, where no measurements are available. The excitation energy is obtained within 
our model from an explicit consideration of the effects of the nuclear potential (Eq. (^) and 
the Coulomb energy (Eq. (P)), i.e. from corrections which are applied to the 4- momenta of 
the final state hadrons leaving the spectator nucleus. We modify the energies of these hadrons 
by the potential barrier and rescale the 3-momenta correspondingly. It is assumed that these 
corrections have to be applied to nucleons wounded in primary and secondary interactions 
and to those hadrons only, which are formed inside the spectator nucleus corresponding to 
the sampled formation time. Among these particles we find apart from the nucleons a small 
fraction of other baryons, which are assumed to move in a nucleon potential and mesons 
to which we apply an effective meson potential of 0.002 GeV. Due to energy-momentum 
conservation these corrections lead to a recoil momentum and, therefore, to an excitation of 
the residual nucleus. In addition, there is a further contribution to the recoil momentum of the 
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residual nucleus arising from potential corrections applied to the momentum of the projectile 
hadron entering the nuclear potential and from cascade nucleons with kinetic energies below 
the nuclear potential which are therefore not able to escape the spectator nucleus. 

In Fig.|l|a we show the dependence of the average excitation energies of the target residual 
nuclei on the momenta of the projectile in the laboratory. The decrease of the excitation 
energy for momenta below about 50 GeV/c is mainly due to the breakdown of the Glauber 
cascade^ as it can be clearly seen for Au and Pb targets. This is not the case for light nuclei 
where even at high energies only up to 2-3 target nucleons are involved on average in the 
primary interaction. This threshold behaviour will be discussed with respect to experimental 
information on grey and heavy particles in more detail further below. At high energies the 
average excitation energies are almost independent on the projectile momentum as one would 
expect from limiting fragmentation [Q. The average excitation energies per nucleon of the 
residual target nucleus are given in Fig.|l|b. In difference to Fig.|l]a the pLab-dependence is 
similar for all target nuclei apart from a constant shift towards higher excitation energies 
per nucleon for light nuclei. This is due to a smaller ratio of wounded nucleons to all target 
nucleons for heavy nuclei as compared to light nuclei. 

In Fig.0 we show the average excitation energy of the residual target nucleus depending 
on the mass number of the target for proton-nucleus interactions at 300 GeV/c. The different 
symbols correspond to several numbers of nucleons lost by the target in primary and secondary 
interactions, i.e. AA is defined by AA = A — A^-es- The excitation energy is strongly correlated 
to the number of removed nucleons. The more nucleons one removes from the target, the more 
energy is deposited into the spectator nucleus. For a fixed number AA the excitation energy 
is increasing with the mass number of the target. The reason for this is, that in heavy targets 
we need more cascading to remove a given number of nucleons than in light targets. As 
an example, the distribution of excitation energies and of excitation energies per nucleon 
of the target prefragment are shown for proton-gold interactions at 300 GeV/c in Fig.^,b. 
In addition to the distributions obtained taking all prefragments into account (labelled "all 
^res") we give distributions which correspond to several mass ranges, in particular to several 
lower cuts in the prefragment mass. Again, the more nucleons are involved in the primary 
interaction and the intranuclear cascade the higher are the mean excitation energies and, 
therefore, the broader are the excitation energy distributions. 

3 Evaporation/Fragmentation 

At the end of the intranuclear cascade the residual nucleus is supposed to be left in an 
equilibrium state, in which the excitation energy U is shared by a large number of nucleons. 
Such an equilibrated compound nucleus is supposed to be characterized by its mass, charge, 
and excitation energy with no further memory of the steps which led to its formation. The 
excitation energy can be higher than the separation energy, thus nucleons and light fragments 

^Note, that the Glauber cascade as obtained with Glauber's formalism is biased by sampling the actual 
chain systems [Q. In order to ensure that the chain masses Af^hain ~ s^p^t exceed the masses of the lowest- 
mass hadronic states with the corresponding quantum numbers lower x-cuts are imposed for all parton systems. 
Therefore, at low energies these x-cuts may reduce the number of sea quark containing chains. 
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(a,d,^H, ^He) can still be emitted: they constitute the low-energy (and most abundant) part 
of the emitted particles in the rest system of the residual nucleus, having an average energy 
of few MeV. The emission process can be well described as an evaporation from a hot system. 
The treatment starts from the formula of Weisskopf ||25| , that is an application of the detailed 



balance principle. The evaporation probability for a particle of type j , mass nij, spin Sj ■ h, 
and kinetic energy E is given by 

PAE)dE = Q^l±^„,jJMEdE (11) 

7r2r Pi{Ui) 

where p's are the nuclear level densities {pf{Uf) for the final nucleus, Pi{Ui) for the initial 
one), Ui = U is the excitation energy of the evaporating nucleus, Uf = U — E — Qj that of the 
final one, Qj is the reaction Q for emitting a particle of type j from the original compound 
nucleus, and cxinv is the cross section for the inverse process. 

Eq. (pTf ) must be implemented with a suitable form for the nuclear level density and the 
inverse cross sections. Many recipes have been suggested for both. In the original work of 
Dostrovsky p6|, p{U) ~ C exp (2-\/ all), with a = A/8 has been used for the level density 
dependence on the excitation energy U. This has led to a simple form for the evaporation 
probability: 

P,{E)dE = + EdE. (12) 

In the same work, the inverse cross sections have been parametrized in a very simple way, 
so that expression (|T2|) can be analytically integrated and used for MC sampling. The same 
formulation is used in this work with, however, a different choice of a as it will be discussed 
later. 

The total width for neutron emission can be found by integrating Eq. (jl^) between zero 
and the maximum possible ejectile energy {U — Qj) 

r,Jl^^f"-''\„„(E)eiE,E. (13) 

Ti'^Ti Jo Pi 

The same applies to charged particles, where the integration actually goes from some effective 
Coulomb barrier where (Tinv drops to zero, up to the maximum energy. 

The evaporative process is in competition with another equilibrium process, that is fis- 
p7| . For the fission probability, a statistical method can be used p5| , |2^ : obtaining for 



sion 



the total fission width 

1 1 r(U~Bp) 



/ pf(U -Bf- E)dE (14) 

Jo 



2np,{U)Jo 

where Bp is the fission barrier, and Pf{Uf) ~ C exp {2^/aFUF) , the level density of the fis- 
sioning nucleus at the saddle point, where the excitation energy Up is given by the initial one 
minus the fission barrier. 

We follow the prescriptions of Atchison |2^ to calculate the quantities entering Eq. (p!^. 



except, again, for the level density parameter ap. 
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In both Pf{U) and p{U) we use the so-called backshifted level density, using U — A rather 
than U, where A is the pairing energy. Moreover, d = a/ A, and dp = a-p/A are found 
to be all but constant parameters: they possess a dependence on A and Z, due to shell 
and deformation effects, and a dependence on excitation energy. Both effects have been 
experimentally observed, and have been subject of many phenomenological and theoretical 



31, 32, 33, 



Here the and Z dependence of Ref. [30] is used 



mvestigations (see 

and complemented with the energy dependence prescription of Ignatyuk 



A-[d-f{U)+d-{l-f{U))] 

ao + 9.17x 10-'' ■ [SziZ) + S^{N)] 



(15) 



fiU) 



0.154-6.3 X IQ-^ ■ A 

I — g-0.054-(;7-A) 

0.054 ■ (f/ - A) 



where according to |^ 



ao IS given by 0.142 MeV~^ and 0.12 MeV"^ for undeformed and 
deformed nuclei respectively, and Sz{Z) and S'n(A^) are the shell correction terms for protons 
and neutrons. The unit of energy used throughout Eq. (|15D is MeV. 

The level density at the saddle point pp is different from that of the nucleus in its ground 
state. From comparison to experimental data, it turns out that ap is greater than the a used 
for evaporation of about 10% at low excitation energies, and the two a's become equal at large 
excitation energies. We use ap ~ 1.08a, with a smooth A dependence. After fission occurs, 
the two fragments are treated like independent residual nuclei with their own excitation and 
can possibly emit further particles. 

For light nuclei, the statistical assumptions and the sequential emission scheme underlying 
the classical evaporation models become less and less applicable, because: 

• Already moderate excitation energies can represent a substantial fraction of the (total) 
binding energy of such nuclei. 

• The level structure of such nuclei is usually highly specific and anyway level spacings 
can be comparable with the excitation energy. 

• The "evaporation" of light fragments other than p or n becomes meaningless, since the 
mass of the "evaporated" fragment can be comparable or even larger than the mass of 
the residual nucleus. 



Therefore other deexcitation mechanisms are more suitable for these light residual nuclei. 
The one adopted for this calculations is the so called Fermi Break-up model 



, ^Tj, where 

the excited nucleus is supposed to disassemble just in one step into two or more fragments, 
with branching given by plain phase space considerations. In particular, the probability 
for disassembling a nucleus of neutrons, Z protons, and U excitation energy (total mass 
M* = U + Ma,z) into n fragments {n > 2) of the same total charge and baryon number, is 
given by: 



br 



n-1 



3/2 



{2tt 



|3(n- 



r 



n 



l)/2 

rp3n/2-5/2 



(16) 
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where the spin factor g, and the permutation factor G are given by {rij is the number of 
identical particles of jth kind) 



n k 



^=n(2^^ + l), G=l[n,l 



(17) 



and is the total kinetic energy of all fragments at the moment of break-up. Vbr is a volume 
of the order of the initial residual nucleus volume. Therefore, the final state are conveniently 
selected by means of a MC procedure, by evaluating such an expression for all possible com- 
binations of fragments energetically allowed and making a random selection. We considered 
all combinations formed by up to six fragments, unless the residual "nucleus" is composed by 
A like particles (p or n), in which case it is disintegrated into A fragments according to phase 
space. All particle stable states with A < 16 have been included, plus the particle unstable 
levels with sizeable 7 decay branching ratios. Also a few known particle unstable isotopes, 
like ^Be, have been included and, if produced, are let to decay according to the experimental 
branching. Once the final state configuration has been selected, the kinematical quantities of 
each fragment are chosen according to n-body phase space distribution. Such a selection must 
be performed taking care to subtract from the available energy the Coulomb repulsion of all 
charged particles: the Coulomb energy is then added back to the charged particles alone, to 
simulate properly the effect of the Coulomb repulsion. In practice E'kin at disassembling will 
be given by: 



where it must be recalled that the emitted fragments can be in an excited state. The total 
Coulomb barrier -Bcoui of the selected configuration is distributed to charged particles after 
disassembling, in their own cm. system. 

According to the picture of the compound nucleus like an equilibrated system determined 
only by its mass, charge and excitation energy, with no memory of previous steps of the 
interaction, Fermi Break-up is activated in the model every time the current compound nucleus 
has mass number A < 17, including possible light fission fragments. The fragmentation of 
higher mass compound nuclei is not yet included in the model. This process, although its 
cross section is quite small, is important when considering the distribution of residual nuclei, 
because it can produce isotopes very far both from the target mass and from the fission 
product distribution. 

4 The production of grey and black particles and resid- 
ual nuclei in high energy collisions 

4.1 Grey and black particles and correlations 

The intranuclear cascade of low energy secondaries and the evaporation of nucleons and light 
fragments mainly contribute to the production of hadrons and light fragments with a velocity 
less than about 0.7c in the rest frame of the target nucleus. They are frequently called "target 




(18) 
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associated" or "slow" particles. Most of the experimental information on slow particle pro- 
duction presently comes from experiments using nuclear emulsions as targets. The emulsions 
usually consist of a component of light nuclei (H,C,N,0) and a component of heavy nuclei 
(Ag,Br). The appearance of slow particles in these experiments has led to their subdivision 
into "grey" and "black" particles. The exact definition usually differs slightly between differ- 
ent experiments. As shown by many authors this subdivision has not only an experimental 
meaning but also subdivides slow particle production into a region which can be understood 
by intranuclear cascade processes (grey particle production) and a region of black particle 
production based on nuclear evaporation processes. If not explicitly stated, throughout this 
paper we apply the following definitions: grey particles are assumed to be singly charged par- 
ticles with a Lorentz-/? value between 0.23 and 0.7 and black particles are singly and multiply 
charged particles with j3 < 0.23. This is in agreement with definitions usually assumed in 
experiments |3^. Furthermore, within our calculation we use the emulsion-composition 



from i.e. an emulsion consisting of 28.8% of light nuclei and of 71.4% of heavy nuclei. 

In Tab.|l] we give the average multiplicities of grey, black, and heavy (=grey+black) par- 
ticles in proton-emulsion interactions as obtained with our MC model DTUNUC 2.0 together 
with experimental results for different momenta of the projectile proton. In addition this 
is shown in Fig.§ together with the corresponding shower particle {(3 > 0.7) multiplicities. 
Whereas the average number of shower particles is increasing throughout the whole energy 
range we get an increasing multiplicity of grey and black particles up to about 40 GeV/c 
which turns into an almost constant behaviour for higher energies. Within our model, this 
constant behaviour is due to limiting fragmentation in each hadron-nucleon interaction [p^ 
together with a constant formation zone intranuclear cascade and inelastic hadron-nucleus 
cross sections depending only weakly on the projectile energy. The position in energy of the 
threshold region, i.e. the region of increasing heavy particle multiplicities, is governed by the 
nuclear geometry and the nuclear potential, which both are treated in a very rough manner, 
and by the way the Glauber cascade is biased by lower cuts applied to chain masses. Further 
shifts of this region or changes of the slope in the threshold region can be obtained by varying 
the distributions from which the x-values of the sea-partons are sampled^. The experimental 
data are taken from a compilation of data by Fredriksson et al. PD|, i.e. they were obtained 



in different experiments. The definitions of "grey" and "black" may therefore slightly vary 
between them. This fact may also partly account for the fluctuations within the data for 
grey and black particle multiplicities in the high energy region. For momenta above about 
20 GeV/c the model agrees well with measured multiplicities, whereas at low energies our re- 
sults seem to depart from the experiments. However, the different definitions used for "grey" 
and "black" cannot completely explain the big differences within the experimental results at 
low energies and any clear experimental information on the threshold behaviour is missing. 
From this comparison it is therefore difficult to draw conclusions about the applicability of our 
model of slow particle production to energies below 10-20 GeV and to modifications within 
the model which could be necessary to reproduce the observed threshold behaviour. 

In order to investigate how the two emulsion components contribute to the energy depen- 
dence of the average multiplicities discussed so far we show in Tabs.Q and ^ and, together with 

* Within DTUNUC they are sampled from a l/a;-distribution H. 
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shower particle multiplicities, in Fig.^ the mean grey, black, and heavy particle multiplicities 
for the light and heavy component separately in comparison to experimental data [Q. In 
order to illustrate the uncertainties within the experimental data we give for some energies 
several multiplicity values, which were measured in different experiments. As it is clearly 
shown in Fig.|^a in interactions of protons with light nuclei even at high energies only up 
to three target nucleons are interacting with the projectile, i.e. we are dealing with a very 
limited Glauber cascade and, therefore, do not observe the typical breakdown of the cascade 
at low energies which would manifest itself in decreasing grey and black particle multiplici- 
ties. In agreement with the measurements we obtain constant mean grey and black particle 
multiplicities in the whole energy range. It seems, that the model overestimates the black 
particle multiplicity which could be due to the absence of the treatment of the nuclear skin in 
the nuclear potential, i.e. by underestimating the low part of the excitation energy distribu- 
tion. However we must note, that the experiments were usually classifying the target nuclei 
as belonging to one of the components by the number of produced heavy particles, which 
implies further uncertainties. An average value of one grey particle per interaction agrees 
well with the experimental results. The model reproduces the measured multiplicities of slow 
particles in interactions of protons with nuclei of the heavy component (Tab.^ Fig.^D) down 
to a proton momentum of about 20 GeV/c. 

In Figs.|^ and |^ we present the grey and black particle multiplicity distributions normalized 
to unity for proton-emulsion interactions at 200 GeV (a) and S~-emulsion interactions at 
350 GeV (b) together with data |Q As the comparisons show, our model is able to 
reproduce the data on slow particle multiplicities very well. The grey particle multiplicity 
distribution for proton projectiles (Fig.^) slightly underestimates the measured distribution 
at high multiplicities which, however, might be not very conclusive since the uncertainties 
within the experimental data are rather big in this region. The hump in the calculated black 
particle distribution for proton projectiles at Nt, ^ 4 (Fig.^a) is due to the evaporation of 
charged particles from light emulsion nuclei and seems to be less pronounced in the measured 
distribution. This is not the case in the S~-emulsion data (Fig.|^), where, on the other hand, 
the uncertainties are higher than within the proton-emulsion data. 

The reasonable description of slow particle multiplicity distributions implies that the model 
should be able to reproduce measured correlations between grey, black, and shower particle 
multiplicities. In Fig.|^ we compare correlations between grey and shower particle multiplici- 
ties (a,b) and between black and shower particle multiplicities (c,d) and in Fig.^,b between 
grey and black particle multiplicities with data of the Alma-Ata-Leningrad-Moscow-Tashkent 
CoUab. |39| on proton-emulsion interactions at 200 GeV and in Fig.^ in addition to data of 



the KLM-Collab. |^2[. In Ref. the errorbars are obviously only given for selected data 
points. Apart from the correlation between grey and black particles, where we obtain slightly 
more black particles for a fixed number of grey particles than seen in the experiments, our 
calculations are in good agreement with the data within their uncertainties. 

A detailed experimental study of slow particle production in interactions of protons, pions, 
and kaons with different target nuclei at energies varying between 50 GeV and 150 GeV was 
presented in Here, grey particles are defined as charged particles having a velocity 

between 0.3c and 0.7c. In Tab.^ we compare our results on mean grey particle multiplicities 
to these data. Again, the agreement is satisfactorily. 
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The dependence of the mean grey, heavy, and shower particle multiphcities on the mass 
number of the target nucleus in proton-nucleus collisions at 300 GeV/c was subject to further 
comparisons. The results are given in Tab.^ and Fig.|l^ together with data taken from Ref. . 



Our model is reproducing the increase of the heavy particle multiplicity with the target mass 
number. As mentioned above the data point for C,N,0 was obtained in emulsion experiments 
in which the identities of the target nuclei were deduced from the heavy particle multiplicities. 
Therefore this data point has to be taken with care. 

Finally, we compare grey particle multiplicity distributions in interactions of protons and 
pions with different target nuclei at 200 GeV/c with recent data of the WA80-Collab. [Q. In 



agreement with the experiment grey particles are defined as singly charged particles with a 
kinetic energy between 30 MeV and 400 MeV emitted in the target rapidity region (—1.7 < 
1] < 1.3). The result of the comparison is shown in Fig.|Tl]. All calculated distributions are 
normalized to the Glauber cross sections of the corresponding interactions. For the two light 
nuclei (C,A1) our distributions are consistently broader than the experimental distributions, 
whereas for heavier targets we agree well in shape and absolute normalization with the data. 

4.2 Residual nuclei and high energy fission 

After evaporation most of the residual nuclei have lost up to one-third of their nucleons 
depending on their mass At, on the kind and energy of the projectile, and on the interaction 
characteristics (impact parameter etc.). They may be considered as heavy fragments produced 
in a spallation or a deep spallation process. In addition, the high energy fission model and 
the Fermi Break-up model which were introduced in Sect.^ modify the mass spectrum of the 



nuclear prefragments furthermore. In Fig.|l^ we show the isobaric mass yields of fragments 



in interactions of silver nuclei with 11.5 GeV (a) and 300 GeV (b) protons together with 



data 0. Since multifragmentation is not included in our MC model we get - apart from 
light fragments {A < 4) which were evaporated from the residual nucleus - almost no fragments 
with masses below Ares ~ 40. In the spallation region (50 < Aj-^s < 100) our calculation agrees 
within a factor of two with the measured mass yields, which is satisfactorily in view of our 
simplified approach and taking into account the fact that multifragmentation would lead to a 
further decrease of the cross section. The rising yields of fragments close to the target mass 
{At — 5 < Aj-cs < At) are not described within our model. This is due to the fact that such 
processes like quasi-elastic scattering are not treated within our model and it might be due 
to our rough treatment of the nuclear potential, i.e. we probably underestimate the low part 
of the excitation energy distribution by neglecting the nuclear skin effects. As experimental 
results on isotope production show, fragment production cross sections remain about constant 
for projectile energies above 10 GeV (see Fig.[l^ and ^). This fact suggests that the 
regime of constant slow particle production may already be reached at an energy of about 
10 GeV. In contrast, within our model the threshold above which slow particle production 
does not change significantly is at about 20-30 GeV (cf. discussion in Sect.^!l|). This fact 
explains the different shape of the calculated mass yields at -^Lab = 11-5 GeV (Fig.|T^) as 
compared to Ej^ah = 300 GeV (Fig.|l2|b) . However, in order to draw further conclusions on the 
threshold region which are based on fragment production cross sections, it would be necessary 
to describe all aspects of the fragmentation process (such as multifragmentation) which is 
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beyond the scope of this work. 

In Fig.0 we compare the charge yield obtained in interactions of 10.6 GeV protons on 
^^^Au nuclei (a) and the isobaric mass yield obtained in interactions of 800 GeV protons 
on ^^^Au nuclei (b) to data [5^]. Since high energy fission significantly modifies the 



fragment production cross sections we show both, the mass yields of the residual nuclei after 
the evaporation-step without high energy fission (crosses) and mass yields obtained taking 
high energy fission into account (diamonds). Within the limitations of our models we are 
able to reproduce the measured yields very well, especially the yields at -^Lab = 800 GeV 
(Fig.|l3|b) where we agree with the data in the mass range 60 < Ares < 190 within their 
uncertainties. The calculated yields at Ares = 2,3,4 represent light fragments evaporated 
from the prefragments. Again, our models do not cover the multifragmentation region and 
the mass region very close to the target mass. 

In order to investigate isotope-production in more detail we compare independent mass 
yield distributions from interactions of 800 GeV protons with ^^''Au with data |ED[ in Fig.in. 



There we plot the cross sections for the production of certain isotopes with masses Ares find 
charge Z versus the difference of their charge and the most probable charge Z^p for three 
intervals of Ares- Corresponding to 0] is defined as 

Z^^{A) = aA^ + bA + c (19) 

with a = -0.382- 10"^ b = 0.483, c = for 82 < Ares < 89, and c = 0.231 for 122 < Ares < 129. 
In the highest mass range 166 < Ares ^ 176 c had to be modified by 1.0 in order to compare the 
shape of the distributions, i.e. c = 1.254. For each interval we calculate the independent yields 
for three different Ares values. We are able to reproduce the measured charge distributions 
which have the typical gaussian shape. 

The average recoil momenta of the fragments in proton-^^^Au interactions at 800 GeV as 
a function of the mass loss AA = At — Ares are shown together with data from different 



experiments |51| in Fig. 15. The momenta of the fragments obtained with our MC model are 



in reasonable agreement with the data. 



5 Summary and conclusions 

We have extended Monte Carlo models based on the Dual Parton Model for high energy 
hadron-nucleus collisions to the calculation of cross sections for residual nuclei production 
and to nuclear evaporation, Fermi Break-up, and high energy fission. 

As it has been demonstrated in a number of past studies 0,]^, ^, Q the models used agree 
quite well with momentum distributions and multiplicities of hadrons produced in high energy 
interactions. Here we find in addition a quite good agreement of the average numbers of grey 
prongs {Ng) and black prongs (Nb) as function of the collision energy and as function of the 
target nucleus with experimental data, which were mostly obtained in emulsion experiments. 

Furthermore, calculated multiplicity distributions of grey and black prongs agree well to 
data. The correlations between the fast shower particles and grey and black prongs as well as 
the correlations between grey and black prongs are often used to analyze the observed events 
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in terms of centrality of the collision or in terms if the impact parameter. Our Monte Carlo 
events show all of these correlations in good agreement with experimental results. 

We find a reasonable agreement of the calculated cross-sections with data for the pro- 
duction of residual nuclei in most of the mass-region below the mass of the original target 
nucleus. Since our model is formulated only in terms of average nuclear properties we can not 
reproduce all the fluctuations, which are due to particular properties of individual nuclei. 

At high energies we find the average numbers of grey and black prongs to become indepen- 
dent from the collision energy. This is a behaviour which can be traced back to the limiting 
fragmentation property of hadron-hadron collisions in the target or projectile rest frame. The 
threshold region, where this high energy behaviour is reached is difficult to predict in a model 
hke ours. As the model has enough freedom to adjust the threshold behaviour to the be- 
haviour of the data, a further tuning of the model parameters might be possible as soon as 
more consistent data become available. 
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Tables 



Table 1: Multiplicities of grey (A^^^), black (A'';,), and heavy {Nh = Ng + Nb) particles in 
interactions of protons with emulsion nuclei are given for different momenta of the incident 
proton. The values as obtained with our model DTUNUC are compared to data from various 



experiments . Within our results we define "grey particles" as particles with a velocity /? = 
v/c between 0.23 and 0.7 and therefore "black particles" as particles with /3 < 0.23. Within 
the experimental data an upper /3-limit for grey particles of 0.7 is usually assumed, whereas 
the /3-cut between "grey" and "black" may slightly differ between different experiments. 



PLab 

(GeV/c) 


(A 

DTUNUC 


Exp. 


(A 

DTUNUC 


^6) 

Exp. 


(A 

DTUNUC 


Exp. 


6.2 


2.0 


3.58±0.11 


3.2 




5.2 


9.25±0.18 


9.0 


2.3 


3.1 ±0.4 


3.6 


4.7 ±0.5 


5.9 




14.9 


2.4 




4.2 




6.6 


8.4 


21.0 


2.6 


2.9 ±0.2 


4.6 


4.6 ±0.2 


7.2 




24.0 


2.5 


3.17±0.1 


4.5 




7.0 


7.7 ±0.2 


50.0 


2.8 


3.07±0.1 


5.0 




7.8 


7.5 ±0.2 


67.0 


2.9 


2.5 ±0.1 


5.2 


4.7 ±0.2 


8.1 




67.0 




2.85±0.09 








7.5 ±0.2 


200.0 


2.9 


2.48±0.08 


5.4 


4.79±0.12 


8.3 




300.0 


2.9 


2.6±0.2 


5.4 


5.4 


8.3 


7.1 ±0.2 


400.0 


2.9 




5.4 




8.3 


8.1 ±0.2 
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Table 2: Grey (Ng), black (A^;,), and heavy {Nh = Ng + Nb) particle multiplicities in inter- 
actions of protons with light emulsion nuclei (C,N,0) are given for different momenta of the 
incident proton. Results of the model are compared to data from various experiments ||40|| . 
For the definition of grey and black particles we refer to the caption of Tab. |1|. 



Phah 


/ A 

DTTTNTTP 

U 1. U IN U 


T \ 
9) 

Exp. 


/ AT 
DTTTNTTP 


\ 

b) 

H/Xp. 


/ A 

TITTTNTTP 
U 1. U In U 


T \ 
h) 

Exp. 


o.U 


0.9 


2.54±0.06 


3.2 




4.1 


8.05±0.1 


o.(J 




U.96±0.07 








2.67±0.14 


9.U 


0.86 


"1 /I 1 n 1 

1.4±0.1 


3.1 




3.96 


4.7±0.1 


zi.U 


1 no 


0. /±0.1 


Z.9 


o_i_n 1 


no 
3.92 


n_i_n 1 
2.9±0.1 


26.0 


1.05 


0.91±0.04 


2.9 




3.95 


2.5±0.1 


50.0 


1.03 


0.91±0.04 


2.8 




3.83 


2.4±0.1 


60.0 


1.07 


0.63±0.07 


2.9 


2.0±0.2 


3.97 


2.6±0.2 


67.0 




0.6±0.1 




1.5±0.1 






69.0 


1.05 


0.84±0.04 


2.8 




3.85 


3.47±0.15 


69.0 




0.91±0.05 








3.65±0.1 


200.0 


1.06 


0.9±0.05 


2.8 


1.8±0.08 


3.86 


2.7±0.11 


200.0 












2.61±0.08 


200.0 












2.75±0.1 


300.0 


1.08 




2.8 




3.88 


2.9±0.21 


400.0 


1.04 


0.67±0.04 


2.8 




3.84 


2.47±0.09 
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Table 3: Grey (A^^^), black (A^";,), and heavy (A^^ = A^^, + A^) particle multiplicities in inter- 
actions of protons with heavy emulsion nuclei (Ag,Br) are given for different momenta of the 
incident proton. Results of the model are compared to data from various experiments ||4CI|| . 
For the definition of grey and black particles we refer to the caption of Tab. 



PLab 

(GeV/c) 


(A 

DTUNUC 


Exp. 


DTUNUC 


Exp. 


(A 

DTUNUC 


Exp. 


9.0 


2.8 


4.1±0.5 


4.0 




6.8 


10.2±0.8 


13.8 


2.9 


6.6±0.6 


4.7 




7.6 


16.0±1.4 


21.0 


3.1 


3.9±0.2 


5.3 


5.9±0.3 


8.4 








3.9±0.2 








9.8±0.3 


24.0 


3.2 


5.1±0.6 


5.4 




8.6 


16.0±1.5 


24.0 




3.96±0.13 








9.5±0.3 


26.0 


3.2 


3.3±0.1 


5.4 




8.6 


11.2±0.15 


50.0 


3.5 


3.86±0.13 


6.0 




9.5 


9.4±0.3 


60.0 


3.6 


3.4±0.2 


6.3 


4.9±0.6 


9.7 


8.3±0.6 


67.0 


3.5 


3.4±0.2 


6.1 


6.2±0.3 


9.5 




67.0 




3.1±0.1 








9.7±0.3 


200.0 


3.8 


3.29±0.1 


6.6 


6.36±0.16 


10.4 


9.66±0.24 


200.0 












9.92±0.17 


300.0 


3.9 




6.8 




10.7 


9.9±0.5 


400.0 


3.9 




6.9 




10.8 


12.4±0.9 


400.0 




3.8±0.1 








9.9±0.2 
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Table 4: Average grey particle multiplicity for proton, pion, and kaon interactions with nuclei 
at different energies. The data are from |^3|. Grey particles are defined as charged particles 
with a velocity v = Pc between 0.3c and 0.7c. 







ELab=50 GeV 


ELab=100 GeV 


ELab=150 GeV 






DTUNUC 


Exp. 


DTUNUC 


Exp. 


DTUNUC 


Exp. 


p 


-C 


0.74 


0.91±0.05 


0.73 


0.82±0.04 


0.77 


0.89±0.04 


p- 


Cu 


2.2 


2.23±0.11 


2.2 


2.26±0.11 


2.3 


2.27±0.11 


p- 


Pb 


4.6 


4.04±0.2 






4.5 


3.75±0.19 




-C 


0.69 


0.85±0.04 


0.69 


0.81±0.04 


0.64 


0.84±0.04 


71+ 


-Cu 


1.96 


1.99±0.1 


1.98 


2.04±0.1 


1.97 


1.99±0.1 


7r+ 


-Pb 


4.0 


3.42±0.17 


3.9 


2.89±0.14 


3.9 


3.31±0.17 


K- 


^-C 


0.65 


0.81±0.04 






0.67 


0.80±0.04 


K+ 


-Cu 


1.82 


1.92±0.1 






1.76 


1.93±0.1 


K+ 


-Pb 


3.62 


3.43±0.17 






3.77 


3.23±0.16 
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Table 5: Dependence of the average multiplicities of grey (A^^) and heavy (Nh) particles on 
the target mass number in proton-nucleus interactions at pLab = 300 GeV/c. Results of the 
model are compared to data from various experiments |40]. "Grey" and "heavy" are defined 
as given in the caption of Tab. |1] 



Nucleus 








DTUNUC 


Exp. 


DTUNUC 


Exp. 


C,N,0 


1.1 




3.9 


2.9±0.21 


Al 


1.4 




5.4 




Cr 


2.1 




7.3 


7.2±0.7 


Emulsion 


2.6 


2.6±0.2 


7.9 


7.1±0.2 


Ag,Br 


3.2 




9.7 


9.9±0.5 


W 


5.0 




13.5 


12.9±1.2 


Au 


5.4 




14.2 
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Figure Captions 



1. Average excitation energies of residual target nuclei in proton-nucleus interactions before 
evaporation are shown for different momenta of the incident proton (a). In (b) the 
average excitation energies are given per nucleon of the residual target nucleus. 

2. Target mass dependence of the average excitation energies for residual target nuclei 
with mass Ares in proton-nucleus interactions before evaporation. AA is the number of 
nucleons lost by the target nucleus in the primary collision {A A = At — Ares)- 

3. In a) the distribution of excitation energies of gold prefragments in proton-gold in- 
teractions at 300 GeV/c is shown for different ranges of prefragment mass ^res- The 
corresponding excitation energy distributions per nucleon of the gold prefragment are 
given in b). 

4. Mean multiplicities of shower, grey, and heavy particles in collisions of protons with 
emulsion nuclei. Data from various experiments |^0| (points) are compared to results of 
the model (lines). 

5. Mean multiplicities of shower, grey, and heavy particles in collisions of protons with 
emulsion nuclei are shown for the component consisting of light nuclei (C,N,0) (a) and 
the heavy component (Ag, Br) (b). Data from various experiments (points) are 
compared to results of the model (lines). 

6. Grey particle multiplicity distributions for interactions of protons (a) and S~-hyperons 



(b) with emulsion nuclei are plotted together with experimental results ^8], ^ 



7. Black particle multiplicity distributions for interactions of protons (a) and S -hyperons 
(b) with emulsion nuclei are plotted together with experimental results I^S], p]] . 

8. The correlations between grey (A^^^) and shower (Ng) particle multiplicities (a,b) and 
black {Nfy) and shower particle multiplicities (c,d) in interactions of protons with emul- 
sion nuclei are compared to experimental results 



9. The correlations between grey {Ng) and black (A";,) particle multiplicities in interactions 
of protons with emulsion nuclei are compared to experimental results ||39|, |42[|. 



10. Target mass number dependence of the average multiplicities of shower, grey, and heavy 
particles in proton-nucleus interactions. Data from various experiments (points) 
are compared to results of the model (lines). 

11. The distributions of grey particle multiplicities in proton-nucleus (a) and pion-nucleus 
(b) interactions as calculated with DTUNUC are compared to experimental results of the 
WA80-Collab. El. 



12. Mass distributions of prefragments produced in proton-silver interactions at 11.5 GeV 
(a) and at 300 GeV (b) as obtained with the model are compared to experimental 



results ^ 
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13. Charge distributions of prefragments produced in proton-gold interactions at 10.6 GeV 
(a) and mass distributions of prefragments produced in proton-gold interactions at 
800 GeV (b) as obtained with the model are compared to experimental results |50[| . 
In addition, the distributions obtained without high-energy fission are shown (crosses). 

14. The yield distributions from interactions of protons with ^^^Au nuclei at 800 GeV are 



shown together with experimental results of Sihver et al. |^0[ for three different intervals 
of the mass of the prefragment. 

15. The total momentum of residual nuclei as a function of the mass loss of the target 
nucleus are compared to experimental results. The experimental data are from different 
experiments and have been taken from Fig. 7 in . 
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p-Nucleus, average excitation energy of the target prefragment 
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p-Emulsion, Ei,ab = 200 GeV 
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P-A, pLab = 200 GeV/c 
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p-Ag, Ei,ab = 300 GeV, mass-yield 
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